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In a major extension of our previous model [^] of condensate growth, we take account of the 
evolution of the occupations of lower trap levels, and of the full Bose-Einstein formula for the 
occupations of higher trap levels. We find good agreement with experiment, especially at higher 
temperatures. We also confirm the picture of the "kinetic" region of evolution, introduced by Kagan 
et al., for the time up to the initiation of the condensate. The behavior after initiation essentially 
follows our original growth equation, but with a substantially increased rate coefficient . 
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Although the hrst Bose condensed atomic vapor was 
produced in a magnetic trap only in 1995 [|[j4|, the ki- 
netics of condensate formation has long been a subject of 
theoretical study There is now intense theoretical 

work on Bose-Einstein condensation, which is excellently 
summarised in j|] . Most theoretical studies of condensate 
growth either have not treated trapping, or have consid- 
ered only traps which are so broad that the behavior of 
the vapor is not essentially different from the untrapped 
situation. Furthermore, they have given only qualitative 
estimates of condensate growth. Our previous paper |^] 
introduced a simplified formula for the growth of a Bose- 
Einstein condensate, in which growth resulted from stim- 
ulated collisions of atoms in a thermal reservoir, where 
one of the atoms enters the lowest trap eigenstate, whose 
occupation thus grows to form a condensate. We thus in- 
cluded the trap eigenfunctions as an essential part of our 
description, and gave the first quantitative formula for 
the growth of a condensate. The growth rate was of the 
order of magnitude of that estimated from experiments 
current at that time. 

This direct stimulated effect must be very important 
once a significant amount of condensate has formed, but 
in the initial stages there will also be a significant number 
of transitions to other low lying trap levels whose popu- 
lations will then also grow. As well as this, there will be 
interactions between the condensate, the atoms in these 
low lying levels, and the atomic vapor from which the 
condensate forms. This paper will extend the descrip- 
tion of the condensate growth to include these factors, 
and will compare the results with experimental data on 
condensate growth obtained at MIT , with which good 
agreement is found. 

As in our previous work, we divide the states in the 
potential into the condensate band, Rq, which consists of 
the energy levels significantly affected by the presence of 
a condensate in the ground state, and the non-condensate 
band, i?NCj which contains all the remaining energy lev- 



els above the condensate band. The division between the 
two bands is taken to be at the value, emax- The situation 
is illustrated in Fig. 1. The picture we shall use assumes 
that i?NC consists of a large "bath" of atomic vapor, in 
thermal equilibrium, whose distribution function is given 
by a time-independent equilibrium Bose-Einstein distri- 
bution {exp[(i? — ^)/kBT] — 1}^^ with positive chemical 
potential /i. The value of Cmax will be assumed to be 
small enough for the majority of atoms to have energies 
higher than emax, so that this part of the bath can be 
treated as being undepleted by the process of condensate 
growth. 
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Fig. 1: a) Representation of the change in the energy spectrum 
due to the growth of the condensate, b) Occupations of the levels 
considered - leftmost is the condensate level, followed by several 
discrete energy bands, with the constantly occupied R^q at higher 
energies. 

In the following we will use the notation of m as the 
mass of an atom with s-wave scattering length a, and no 
for the number of condensate atoms. 

The levels in i?c must have time dependent energies 
due to the effects of the interaction with the growing 
condensate. The energy of the ground state is /x(no), 
the chemical potential, for which we use a modified 
Thomas-Fermi approach /i(no) = ain^ -\- v)'^/^ where 
a — {15aujj;UjyUJz''n^^'^fi^/4:V^)'^^^ and ly is a constant 
chosen so that /i(0) = h,{u!x -\- tOy -\- ujz)l2. As ^{uq) 
rises with an increase in no , the energies of higher energy 
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levels must also rise. The exact nature of this rise does 
not affect the results greatly. An approximate treatment 
arises by leaving the levels with energies above 2/i(no) 
unchanged (including all levels in i?Nc)j and compress- 
ing the spacing of the levels under 2fi{no), so that the 
energies are given by 



?(2^(no)-e°„)(2MK) 



(1) 



where are the non-interacting harmonic potential en- 
ergy levels, and 9{x) is the step function. Note that 
Cq = /i(0), and thus eg — /i(no). The levels used in this 
model are represented graphically in Fig. 1. To simplify 
the equations we also group the levels in narrow bands 
of mean energy e^, gk levels per group, and Uk particles 
per group. This corresponds to the ergodic assumption 
used in 1,0,0 

Equations of motion then follow from Quantum Ki- 
netic theory [pl and full derivation will be given else- 
where, but the result is essentially equivalent to modify- 
ing the Quantum Boltzmann (QBE) equation as follows. 

(i) We use the QBE in an approximated ergodic form, 
where Cmin = min(e™, e„, e^, e,) [|| 



5/(e„) 



[/(ep)/(e,)(l + /(e™))(l + /(e„)) 
-/(e„0/(e„)(l + /(ep))(l + /(e,))l, (2) 



where Uk = gkf{&k) is the number of particles with en- 
ergy Cfc, and T = {Sma'^uj'^ /nh) x t. 

(ii) We use the modified energy levels as given above, 
but otherwise do not change the QBE. 

(iii) We sum out over all levels in i?NC which is as- 
sumed time independent. 

(iv) We omit any scattering between atoms which both 
have energies less than emax, which is reasonable if the 
number of atoms in the bath is almost 100% of the total 
number of atoms. 
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Fig. 2: The transitions being considered: Left — scattering; Right — 
Growth. 

As a result of this procedure there are two different 
kinds of dynamics, scattering and growth, as illustrated 



in Fig. 2. The evolution for the population of the mth 
level in Rq is then 



dr 



dr 
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m I growth ~r ''^'mlscatt 



(3) 



where the individual terms are as follows. 
Scattering: A collision between an atom initially in an 
energy level below Cmax and a bath atom transfers the 
first atom to another energy level below Cmax- This is 
described by 



— nfc(ffm + n„)e ''"^"-"^'"''^ -n„,{gk + nk) 

k<m ^"^ 

- V — Ufe(.g„ + Tim) - n^igk + nfe)e"''"'""/'''«^ 



(4) 



where LOkm ~ Sk — em- The value of r(T) = 
g-fim/fcisT depends on the energy spectrum. 
We use the value for an isotropic 3-dimensional har- 
monic oscillator with frequency uj = {ujxtOytOzY^^ ■ Thus 
e„ = (n + 'i/2)TiLo, so that we find 



r(r) 



(5) 



However this value is not critical; similar results are ob- 
tained for any r(r) greater than about 10% of (H). 
Growth: A collision between a pair of atoms initially 
in the bath of atomic vapor results in one of the atoms 
having a final energy less than Cmax- This gives 



I growth 

where 



2[inm + l)W+{no)-nmW-{no)] (6) 



dei 



dey 



des.f{ei)f{e2) 



X (1 + /(e3))<5(ei -I- 62 - 63 - e™). 



(7) 



W-{n^) - - 



OC /'OC 

de-) 



dei J de2 J ^63(1 -|- /(61)) 

X (1 + /(62))/(63)<5(6i + 62 - 63 - e^). (8) 



Because we assume the non-condensate band is in ther- 
mal equilibrium, with temperature T and chemical po- 
tential fj,, we have 

W^^(no) = cxp (%^) W+ino) (9) 

The value of (henceforth differs significantly 
from that used previously [|| in that the integrals in (|^) 
are only performed over energy levels higher than 6inax 
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(previously over all levels), and the Bose-Einstein distri- 
bution function is used over all the range of integration 
(previously approximated by the Maxwell-Boltzmann 
distribution). With these changes we find that 

+z2^[zz(no)r[a>(z,l,r+l)]n, (10) 

r=l ■' 



where 



z(no) 







(11) 



The function $ is the Lerch transcendent , defined by 
$(x, s, a) = J2'k'=o ^'^ I + ■ These changes result in a 
significant correction, increasing by about a factor 
greater than three (dependent on T, and the trap pa- 
rameters) compared to that used previously in Q, and 
producing correspondingly faster growth. 

By making a further approximation, that W^^{no) « 
W^{no), the calculations required are significantly sim- 
plified. We can do this because the iy+(no) are an av- 
erage over all the levels contained in the mth group, and 
hence are expected to be of the same order of magnitude 
as VF+(no). As a validity check, it was found that the 
effect on the condensate growth rate was small when the 
W^(no) were altered by a factor in the range 0.5 — 2. We 
now have, for the growth terms 



J^m I growth = 2VF+(no) I 
?^o I growth = 2W^+(no) 



1 — e 



1 - e ''bt 



(12) 

no + l|. (13) 



The overall evolution of the system can now be found 
from the numerical solutions to the set of equations (|^). 
The parameters used are Wa; = Wj, = 27r x 82.3Hz and 
a;^ = 27r X 18Hz, as in and a — 2.75nm, e,nax ~ 
2.2^(no(eq)) , where no(eq) is the final equilibrium value 
of rig, which is defined by /i(no(eq)) — /i, the chemical 
potential of the vapor. The number of energy bands was 
set as large as possible, but it was required that there 
were at least four levels in the first group above the con- 
densate, in order represent the fact that the levels are 
discrete. 

The improvements to this model over that used in 0] 
speed up the condensate growth by up to an order of 
magnitude, depending on the exact parameters, as antic- 
ipated in 0. The major cause of this is the correction to 
arising from the use of the correct Bose-Einstein dis- 
tribution. The inclusion of the scattering terms does not 
change the overall rate of growth substantially (which 
is dominated by the bosonic stimulated growth), but 
does speed up the initial period of growth (dominated 



by spontaneous growth) and shortens the time before the 
stimulated term becomes significant. This gives a much 
sharper initial growth curve as compared to the smoother 
S-bend curves of In Fig. 3 we present an example of 
the results obtained for the growth of all the bands in 
i?c , in which a number of features may be seen. 

(i) The effect of the initial conditions used can be seen 
from the front corner. In this example the initial pop- 
ulations for all but the top ten bands were set at zero, 
whilst the top ten bands had populations determined by 
the same Bose-Einstein distribution function as for i?NC • 
The figure shows that this initial condition is rapidly 
smoothed out by scattering processes. Different initial 
conditions merely generate a small change in the initia- 
tion time defined below. 
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Fig. 3: f{em) vs. energy Cm as a function of time. Note that the 
lines almost parallel to the time axis are not lines of constant en- 
ergy, but rather lines of /(e,n) for constant level number whose en- 
ergies change with condensate growth. The solid curve in the plane 
at the top of the axes represents the curve log]^Q[eo] = logjQ [/i{rto)] 
as a function of time. 

(ii) The initiation time (here 60ms) is defined by the 
critical line C-C. Up to this point the population of the 
condensate level is relatively small. The behavior up to 
the initiation time is similar to that found by Svistunov 
|]l0| and Kagan et al. Q . In particular the populations of 
the levels increase to approach a limiting dependence on 
energy of f{E) cx E~^-^^ on the critical line C — C, which 
is in good agreement with their prediction of E~^/^ |l^] 
for the case of a harmonic trap. 

(iii) After the initiation time the condensate grows 
enormously. However, the occupations of the other trap 
levels actually decrease quite rapidly to their equilibrium 
values, and then remain nearly constant while the con- 
densate continues to grow - by a factor of about 10 in this 
case. At the same time the energy spectrum of the trap 
levels in Rq changes according to ([^). This accounts for 
the small variation in occupations of these levels which 
is still evident in Fig. 3. 

Fig. 4 compares with the experimental data of ref. 
for two different temperatures. The MIT group fit- 
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ted their data to an uncorrected growth equation no = 
7no[l — (no/no(eq))^/^], and reported only values of the 
parameter 7. To represent that the experimental curve 
is a fitted curve rather than the raw data, the curve has 
been plotted as a broad band. The MIT group used the 
initial population for their curves as free parameters in 
their fit. We have set the initial populations for the MIT 
curves so as to give the best agreement with the initiation 
time of our growth. 
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Fig. 4: Comparison of tlieoretical growtli curves (black) with ex- 
perimentally fitted curves (gray) from jsj. a) Theory T = 830nK, 
no(eq) = 7.6 X lO"; Experiment T = 810 - 890nK, no(eq) = 
7.5 - 7.85 X 106; b) Theory T = 590nK, no(eq) = 2.3 X lO"; 
Experiment T = 580 - 610nK, no(eq) = 2.1 - 2.5 X 10*^. 

In the T ~ 830nK case, Fig. 4a, the growth speed 
predicted agrees with that experimentally found. The 
T = 590nK case in Fig. 4b shows a theoretical growth 
rate some three times slower than is measured. The dis- 
crepancy between theory and experiment at lower tem- 
peratures is hard to evaluate using the data in the form 
presented in which do not allow for direct compar- 
ison between the actual projected spatial distributions 
as given by phase-contrast microscopy, and theoretical 
spatial distributions. These are not difficult to calcu- 
late from our many-level growth curves — the methodol- 
ogy will be published elsewhere. The MIT method fits 
to a zero chemical potential vapor plus a non-zero chem- 
ical potential condensate - a reasonable estimate in the 
absence of any theoretical description of the spatial dis- 
tribution of the vapor. But a detailed description might 
give quite different results for temperature, and for con- 
densate and vapor numbers. 

In summary, we have given a description of condensate 
growth which covers the full range of behaviors, both be- 
fore and after initiation of the condensate. It is quantita- 
tive and agrees quite well with experiment. The behavior 
before initiation is essentially as given by the solution of 
the quantum Boltzmann equation, and agrees with com- 
putations based on this equation But we are 
able to give a value for the initiation time which appears 
to be consistent with experiment. After the initiation of 
the condensate, the occupations of the non-condensate 
levels are clamped by their fast relaxation time to quasi- 



equilibrium values, which change with the rise in trap 
levels induced by the very much slower growth of the 
condensate to its final occupation. 

We believe the future development of the theory for 
this problem will involve mainly refinements of our pic- 
ture, such as including depletion of our fixed "bath" of 
vapor. But for quantitative comparison with experiment, 
it will be necessary to have more extensive data on spatial 
distributions at a range of temperatures. 
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